n = 51;
A = readmatrix("test1-pol.csv"); Z1pol = A(:,1:n)';
A = readmatrix("test2-pol.csv"); Z2pol = A(:,1:n)';
A = readmatrix("test3-pol.csv"); Z3pol = A(:,1:n)';

[X,Y] = meshgrid(1:0.02:2,1:0.02:2);
subplot(1, 3, 1);
surf(X,Y,Z1pol); xlabel('X'); ylabel('Y'); zlabel('Z'); grid on;
view([15, 10])
subplot(1, 3, 2);
surf(X,Y,Z2pol); xlabel('X'); ylabel('Y'); zlabel('Z'); grid on;
view([15, 10])
subplot(1, 3, 3);
surf(X,Y,Z3pol); xlabel('X'); ylabel('Y'); zlabel('Z'); grid on;
view([15, 10])
